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In this paper, we theoretically show that a broadband resonant enhancement of emission may 
occur for infrared sources located in a polaritonic wire medium. The reason of this enhancement is 
overlapping of two topological transitions of the wave dispersion in the medium. The first topological 
transition has been recently revealed as an effect inherent to polaritonic wire media. The second 
one was theoretically uncovered in another material. In this work we reveal it for wire media and 
prove the possibility to combine both these transitions with the purpose to obtain the broadband 
resonant Purcell factor. We compare the results obtained for two orientations of a subwavelength 
electric dipole embedded into wire medium - that along the optical axis and that perpendicular to 
it - and report on the resonant isotropic radiation enhancement. Also, we reveal the enhancement 
of radiation to free space from a finite sample of the wire medium. 

PACS numbers: 41.20.Jb, 81.05.Xj, 78.67.Uh, 74.25.Gz 


I. INTRODUCTION 


Throughout the last decade, enormous efforts have 
been put into research and development of artificial elec¬ 
tromagnetic media with unusual properties. One of these 
media is the so-called wire medium (WM), defined as an 
optically dense array of parallel wires embedded in a host 
dielectric matrix, as depicted in Fig. |l(a)[ Its unique 
feature is applicability from radio frequencies to those 
of the visible light (see, e.g., Ref. i)- Subwavelength 
imagiM 0 Q , improvement of thermophotovoltaic sys¬ 
tems [J, enhancement of directionality of microwave 
antennas are some examples of such applications. One 
of prominent applications is optical sensing and location 
of fluorescent nano-objects at optical frequencies. Gen¬ 
erally speaking, this possibility is provided because the 
power radiated by the fluorescent emitter located in such 
a medium may increase several times or even by an order 
of magnitude (see, e.g., Ref. 0)- 

More in detail, since wire medium is optically dense, it 
can be homogenized, and based on the effective medium 
model, it is a uniaxially anisotropic material whose per¬ 
mittivity tensor has principal components which (at least 
at some frequencies) have different signs 0. Therefore, 
the sign of the permittivity tensor trace is indefinite, and 
wire medium can be classified as an indefinite material 
according to Ref. 0. This feature affects remarkably 
the isofrequency surface (dispersion surface in the recip¬ 
rocal lattice called k-space) of the extraordinary modes 
propagating in the medium, whereas the isofrequency for 
the ordinary mode keeps similar to that of an isotropic 
medium. The indefinite regime implies an open surface 
close to a hyperboloid 0. Hence, at corresponding fre¬ 
quencies, wire medium refers to the class of so-called hy¬ 
perbolic metamaterial (see, e.g., Ref. 0). The hyperbolic 
topology implies very high spatial frequencies of propa¬ 
gating extraordinary modes and, respectively, very high 


density of electromagnetic states. This effect results in 
the enhancement of the radiated power of a dipole source 
submerged into the medium (see, e.g., Refs. 0, [HI)- 
Increase of the radiated power implies enhancement of 
spontaneous emission and decay rate. The enhancement 
of the decay rate due to an environment of an emitter 
is called Purcell’s effect. Historically, this effect was re¬ 
ferred as the enhanced decay rate of the excited state of 
a quantum emitter due to insertion of this emitter into 
an open cavity However, recently the concept 

of Purcell’s effect was generalized to any resonant scat- 
terer located in the vicinity of the emitter [IM3 , and 
finally - to the case of an arbitrary environment different 
from free space (see, e.g.. Refs. 0[il, [3). The increase 
of the radiated power of a monochromatic point dipole 
due to its environment is called the Purcell factor. Con¬ 
sequently, WM may offer a high Purcell factor in the 
hyperbolic regime. 


High Purcell factor is inherent to any hyperbolic ma¬ 
terial. A periodic stack of strongly submicron bilayers, 
where one layer is a plasmonic material and another is 
a material with positive permittivity, is a popular im¬ 
plementation of a hyperbolic medium operating in the 
optical frequency range (see, e.g., Refs. i0)- There¬ 
fore, this medium also allows a high Purcell factor [i0 . 
Moreover, the Purcell factor for such medium can be reso¬ 
nant, i.e. attains a maximum at a certain frequency. This 
resonance is claimed for the case when a transition hap¬ 
pens in the topology of the isofrequency surface. Namely, 
the resonant maximum of the Purcell factor occurs at a 
frequency where the isofrequency of the stacked hyper¬ 
bolic medium transits from the closed-surface (elliptical) 
regime to open-surface (hyperbolic) regime or vice versa 
2l|. This transition was called topological phase transi¬ 
tion 


In our recent work [23|, we studied the resonat be¬ 
havior of the Purcell factor for a wire medium formed 
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by nanowires prepared from a so-called polaritonic mate¬ 
rial and operating at infrared frequencies. The resonance 
corresponds also to the topological transition. However, 
the topological transition in this medium, unlike that of 
the stacked material, offers an unbounded spatial spec¬ 
trum of radiation to internal sources. At the topological 
transition revealed in Ref. [l^, the isofrequency surface 
is composed of both the hyperbolic and elliptic disper¬ 
sion branches connected with one another as depicted in 
Fig. |l(b)[ Since the spatial spectrum is unbounded over 
the parallel wave vector (3 (vertical axis in the Fig. |l(b)1 , 
a parallel dipole generates more radiation than the per¬ 
pendicular one [ 23 . This result makes our work Ref. [ 2 ^ 
actual in view of the aforementioned application. 

In the present paper, we report that the same polari¬ 
tonic WM as was studied in Ref. [ 2 ^ experiences, in fact, 
several topological transitions. Two of them are impor¬ 
tant since correspond to the resonant enhancement of the 
dipole radiation. Restricting the frequency band of our 
interest, we have not noticed the second useful transition 
in Ref. [ 2 ^. In fact, the two transitions occur at neigh¬ 
boring frequencies. On some feasible conditions, two cor¬ 
responding resonances of the Purcell factor may overlap. 
This overlapping has an important implication for the 
radiation of dipole sources in a polaritonic WM. Both 
values and frequency dependencies of the Purcell factor 
turn out to be very different from those predicted in the 
literature for a stacked hyperbolic material. This differ¬ 
ence is favorable for a WM and makes it, to our opinion, 
more attractive for infrared nanosensing than the stacked 
hyperbolic materials. This is the main message of the 
present work. 

The paper is organized as follows: in Sec.ITTlwe discuss 
these two topological transitions, in Section IlIIl we sum¬ 
marize the results, and finally in Sec. lIVl we conclude the 
paper. 


II. THEORY 


Figure 1(a) schematically illustrates a WM as a regular 
optically dense array of parallel cylindrical wires embed¬ 
ded in a host dielectric material. In our calculations, we 
assume for simplicity that the host material is free space. 
As the figure shows, the optical axis of the medium has 
unit vector Zq which is also the wire axis. In this case, 
homogeneous effective permittivity tensor of the medium 
is given by 
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FIG. 1. (a)-A medium of parallel wires which form a square 
lattice. (b)-Isofrequency of the first useful topological transi¬ 
tion comprises both closed-surface (elliptic) and open-surface 
(two-branch hyperbolic) surfaces. (c)-The second useful 
topological transition occurs as a sharp (in the lossless case) 
substitution of an ellipsoid by a solid hyperboloid. 
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The subscripts T and || mean components, perpendicu¬ 
lar and parallel to the optical axis, respectively. In Eq. [2l 
/v is the volume fraction and £r represents the relative 
dielectric constant of the wires. Furthermore, fcg and 
kp are the free-space and plasma wave numbers, respec¬ 
tively, and j3 is the parallel component of the wave vector. 
The effective-medium model |24| was developed for wires 
performed of a material with negative dielectric response. 
Its applicability is determined by two main conditions - 
that of sufficiently low spatial frequencies qb < tt {q is the 
perpendicular component of the wave vector) and that of 
long wavelengths kob <C 27r. If the radius of the wires is 
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much smaller than the wavelength {koa <C 1) and the 
fraction of wires is rather small (practically /v < 0.3), 
the effective-medium model [23 ]j in accordance to sev¬ 
eral recent studies, is quite accurate. 

Distinguished by directions of the electric and mag¬ 
netic fields, two different modes exist in the WM. The 
first one is transverse electric (TE) or ordinary mode, 
whose electric field is normal to the optical axis. The 
dispersion equation for this mode is standard: 

q^+ii^ = kle^. (3) 


The isofrequency is a simple sphere with the radius 
If the fraction volume is small and the relative dielectric 
constant of the wires is not close to —1, the perpendicular 
component of the effective permittivity is close to unity. 
It implies that the TE modes do not interact with the 
wires and the Purcell effect for a source creating TE- 
polarized radiation is absent. 

The second eigenmode is transverse magnetic (TM) or 
extraordinary mode whose magnetic field is normal to 
the optical axis. Its dispersion equation can be written 
as follows (see, e.g.. Ref. [^1: 
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As mentioned in the Introduction and Eq. 2] determines, 
depending on the sign of ey and e_L (more exactly, the 
real parts of these components), an open surface or a 
closed surface arises in the reciprocal space. Definitely, 
5i(e||) < 0&5R(ex) > 0, or 5ft(e||) > 0&5R(ej_) < 0 gives 
an open surface, and 5i(£||) > 0&5R(£x) > 0 results in a 
closed isofrequency surface. These surfaces can be differ¬ 
ent from, respectively, hyperboloid and ellipsoid due to 
the spatial dispersion (dependence of the effective per¬ 
mittivity on /3). 

In the following subsections, we discuss two different 
types of transitions from one dispersion surface to an¬ 
other one. At the first type of transition, the two-branch 
hyperboloid transits versus growing frequency to an ellip¬ 
soid. This transition in polaritonic WM comprises several 
stages which may be in their turn called separate topo¬ 
logical transitions. At low frequencies the isofrequency 
surface of our WM is a two-branch hyperboloid. Start¬ 
ing from some frequency there appears an small ellipsoid 
which enlarges versus frequency and coexists with the 
hyperboloid. This is the first stage of the topological 
transition studied in Ref. [l^. The second stage cor¬ 
responds to the certain frequency at which both hyper¬ 
bolic and elliptic surfaces connect one another and form 
a solid isofrequency. At the third stage the cut-off of 
the hyperboloid in the reciprocal space increases and it 
disconnects from the ellipsoid. This growth of the cut¬ 
off versus frequency is very fast, and at some frequency 
the hyperboloid disappears from the domain of physi¬ 
cally sound values of the wave vector. This dynamics of 
the isofrequency surface can be treated as a multi-stage 
transition from the hyperbolic to the elliptic dispersion 


of propagating waves. However, in Ref. (23|, we mainly 
concentrated on the stage when the hyperbolic and ellip¬ 
tic surfaces are connected, since it grants the unbounded 
spatial spectrum of eignemodes. This regime grants the 
strongest enhancement to the dipole radiation and we ob¬ 
served the maximum of simulated Purcell factor namely 
at this frequency. 

At the second topological transition, illustrated by 
Fig. 1(c) an ellipsoid is substituted by a solid hyper¬ 
boloid. This hyperboloid describes the dispersion in uni¬ 
axial media whose perpendicular permittivity is negative 
whereas the parallel one is positive. Below we show that 
this is so for polaritonic WM at sufficiently high frequen¬ 
cies. Therefore, both types of topological transition hold 
in the same material - WM of polaritonic nanowires. 

Notice that in Ref. Q, the resonance of the Purcell fac¬ 
tor of a WM performed of £r-negative nanowires was re¬ 
vealed, though not explained. We have calculated isofre¬ 
quencies for WM whose Purcell factor is depicted in 
Fig. 4 of Ref. 0 and found that the resonance of the Pur¬ 
cell factor for a parallel dipole (which occurs at £r = —90, 
kob/jT = a/b = 0.05) corresponds namely to our topolog¬ 
ical transition of the first type. 


A. Topological transition via unbounded spatial 
spectrum 


Though this transition has been studied in Ref. [^ . 
in this subsection, we briefly review this work for the 
convenience of the reader. 

Equation [3 where we adopt for simplicity an approxi¬ 
mation £_l Ri 1, together with Eq.|3] delivers two solutions 
for TM-waves: 
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Here it is denoted: 
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The first solution corresponds to open-surface disper¬ 
sion (e|| < 0), and the second solution - to closed-surface 
dispersion (e|| > 0). The separation of these two surfaces 
in the reciprocal space is determined by kc- If |£r| is very 
large, kc is much smaller than the plasma wave number 
kp. In this scenario, £|| is negative at the corresponding 
frequency and the isofrquency surface is open suffering 
a cut-off over the parallel component of the wave-vector 
axis: \f3\ > ko. The closed surface in the reciprocal space 
is absent since the second solution in Eq.[5]corresponds to 
/3^ < 0. However, if |£r| is a modest value, the parameter 
kc becomes comparable to kp. Now, the second solution 
like the first one also corresponds to /3^ > 0. Therefore, 
the closed surface also appears and both surfaces coexist 
as parts of the complex isofrequency surface. Note that 
the closed surface covers low values of the wave vector: 
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1^1 < fco- Finally, at the frequency where kc = fcp, the 
first (1/31 > ko) and the second (|/3| < fco) solutions touch 
each other at /3 = i/co and the surfaces are connected. 
This connection gives topological phase transition of the 
first type which is the most important stage of this multi¬ 
stage process because the spatial spectrum becomes un¬ 
bounded. Not in any WM the regime with transition of 
kc can approach to fcp is achievable, however, such WM 
exists. 

At the points /3 = ±fco where the open and closed sur¬ 
faces are connected, the parallel component of the effec¬ 
tive permittivity is zero e|| = 0. As mentioned, fcc = fcp, 
therefore, after some easy algebra, we obtain the condi¬ 
tion of this desired regime in the form 



This condition is sufficient for the topological transition 
of the first type in the framework of the effective-medium 
model [23 • 

Suppose that the frequency dependence of the dielec¬ 
tric response of the wires can be modeled as 


Erito) = — 


fco 


Ki 

K2 - jKs ’ 


( 8 ) 


For low-loss materials, we may approximately equate 
right-hand sides of Eq. [7] and Eq. [S] Optical losses in the 
material of nanowires, if sufficiently small, are irrelevant 
for this topological transition, and in this subsection we 
neglect them. Of course, in our work [ 2 ^ we took these 
losses into account and considered their influence. 

In Fig. [21 the normalized perpendicular wave vector 
is plotted versus the normalized parallel wave vector 
for four normalized frequencies fcp6 = 0.9, fcp6 = 0.95, 
kob = 1, and kob = 1.05. In these calculations, we put 
a/b = 0.05, £(“) = 1, Kib = 50, ATzb = 0.6073, and 
K 3 = 0. At kob = 0.9, the isofrequency surface (magenta 
color) is open and close to a hyperboloid. This frequency 
corresponds to = —169.82 and fcc/fcp = 0.8633. At 
higher frequencies |£r| decreases and fcc approaches fcp. 
At kob = 0.95 (the first stage of the first-type topologi¬ 
cal transition) the open and closed branches coexist but 
are separated. Their touching occurs at kob = 1 when 
£ 1 . = 1 — l//v = —126.324 and fcc = fcp. At this fre¬ 
quency the spatial spectrum is unbounded. This is the 
second and most important stage of the first-type topo¬ 
logical transition. At kob > 1, |£i.| further decreases and 
the two surfaces separate again (the third stage). At a 
higher frequency the hyperbolic branch disappears and 
only the elliptic one remains (the fourth stage). 


B. Topological transition from ellipsoid to solid 
hyperboloid 

If we consider an extended frequency range, we may 
notice that the another topological transition occurs in 
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FIG. 2. Isofrequency contours (cross section of the dispersion 
surface) for TM mode. Calculations are done based on the 
following values: a/b = 0.05, Kib — 50, K 2 b = 0.6073, K^b = 
0 and = 1. Magenta color corresponds to normalized 
frequency kob = 0.90, red color - to kob — 0.95, blue - to 
kob = 1, and black ~ to kob = 1.05. 


the same WM. Really, based on Eq. [21 we may write 

1 £r(a;) — Z 
^ ~Zec{uj)-P' 

where the pole P and zero Z are expressed as 


(9) 


l + /v 1 - /v 

l-/v’ l + /v’ 


( 10 ) 


Both the zero and the pole of £_l exist at higher frequen¬ 
cies than that of the first topological transition. As the 
above equation indicates, P < —1 and —1<Z<0. 
If the fraction of nanowires is very small, the pole and 
zero are close to each other on the frequency axis. Defi¬ 
nitely, at the pole (£r = P) and at the zero (£r = Z) the 
sign of the perpendicular component e± changes. Recall, 
that the first topological transition occurs at a frequency 
where Re(£||) may vanish (vanishes when |/3| = fco). The 
zero of Re(£_L) also results in a topological transition. 
This second transition exploits the frequency dispersion 
of £r. The spatial dispersion is here insignificant - this 
zero occurs at the corresponding frequency for any /3. 

Definitely, the hyperbolic isofrequency extends to in¬ 
finity and therefore implies the broader spatial spectrum 
than the elliptic one. Based on this one could expect the 
growing Purcell factor at frequencies above this topo¬ 
logical transition. However, in fact, a local maximum of 
Purcell fcator occurs exactly at this frequency. This max¬ 
imum was explained in Ref. [ 2 ^ for the case of a stacked 
hyperbolic medium using the previously obtained Pur¬ 
cell factor of a uniaxial dielectric. Really, the zero of 
the perpendicular effective permittivity, which occurs at 
the topological transition of the second type (whereas 
the parallel permittivity keeps finite), implies the ulti- 
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FIG. 3. Perpendicular component of the effective permittivity 
(real and imaginary parts - solid and dashed, respectively) 
versus normalized frequency. Red and black curves corre¬ 
spond to ajh = 0.05 and a/6 = 0.08, respectively. Calcula¬ 
tions are done for the case = 1, K\h — 50, K 2 b = 0.6073 
and Ksb — 0.03. 


mate anisotropy of the medium. Consequently, the ra¬ 
diation of any finite source (whatever small) should be 
super-directional. The super-directionality of the dipole 
radiation implies the singularity for the density of pho¬ 
tonic states, and the Purcell factor becomes maximally 
high. If both parallel and perpendicular components of 
the effective-medium permittivity tensor are positive the 
formula for the parallel-dipole Purcell factor describing 
this effect takes form [^ : 


Pf 


g|| 



( 11 ) 


This estimation gives an the infinite value of the Purcell 
factor in the lossless case, because does not take into ac¬ 
count the limit of the effective-medium model for spatial 
frequencies. Definitely, this formula is not applicable for 
lossy hyperbolic media, it only gives an insight of the 
effect of ultimate anisotropy e|| ^ e±. 

Fig.Oshows e_L versus frequency, where Eq.[5]describes 
the frequency dispersion of the material of nanowires. In 
the present calculations, the losses in this material are 
taken into account: parameters of the material are cho¬ 
sen as follows: = 1, Kib = 50, K 2 b = 0.6073 (as 

before), and K^b = 0.03. Optical losses transform the 
pole and the zero of the perpendicular permittivity into 
a Lorentz’s type resonance. The pole is substituted by 
the maximum of |ex| and the zero refers in the lossy case 
to the real part. An increase of the volume fraction of 
nanowires causes the red shift of the resonance frequency 
ex and enhances the magnitude of the resonance. Notice 
that in the whole frequency range of Fig. [31 the parallel 
component of the effective permittivity e|| is positive (at 


spatial frequencies for which the effective-medium model 
is valid). Therefore, the transition from the ellipsoid 
to solid hyperboloid in the case a/b = 0.05 occurs at 
kob ~ 25.4, and in the case a/b = 0.08 - at kob ~ 25.2. 
Unlike our previous subsection, we do not depict here the 
dynamics of the topological transition because it is sim¬ 
ilar to that in the stacked hyperbolic metamaterial (see 
Fig. 3 of HI). 

At frequencies kob ~ 26 {a/b = 0.08) and kob = 25.7 
{a/b = 0.05), the solid hyperboloid is again substituted 
by an ellipsoid. This topological transition is also char¬ 
acterized by 5ft(ex) = 0 and implies the maximum of the 
Purcell factor. However, as we will see below, in practi¬ 
cal lossy media these two maxima do not exist separately. 
The ratio e||/ex may be large enough at the local mini¬ 
mum of ex(w) which occurs at a Lorentz’s resonance of 
the effective medium without crossing the zero. Then 
the Purcell factor also experiences a resonance though 
the dispersion keeps elliptic. This regime (mimicking the 
2 d topological transition by the minimum of ex(t<-')) will 
be also considered in the next section. 

Since both topological transitions in a WM - that of 
the first type and that of the second type - correspond to 
the resonances of the Purcell factor, an interesting prob¬ 
lem arises: to approach these two resonances and, over¬ 
lapping them, in order to obtain a broadband and strong 
Purcell effect where the values of the Purcell factor would 
dramatically exceed those obtained in Refs. HI and HI 
for stacked metamaterials and even the values obtained 
in 0 for perfectly conducting WM. In the following sec¬ 
tion we design the needed WM and theoretically reveal 
the expected effect for its finite-size sample. 


III. OPTIMIZATION OF THE WIRE MEDIUM 
SAMPLE 


We start from the volume fraction of nanowires /v_= 
0.0804. This value was used in our previous study H3 
in order to demonstrate the topological transition of the 
first type. In this case, Eq. |7| gives (assuming that the 
optical losses are sufficiently small) the relative dielectric 
constant of the wires 5ft(e:r) = —11.4. This dielectric re¬ 
sponse is modestly negative. Polaritonic materials like 
lithium tantalate (LiTaOa) can provide such response at 
infrared frequencies. For this range, the chosen volume 
fraction may correspond to a lattice constant b =200 nm 
and wire radius a =32 nm which are much smaller than 
the wavelength. Therefore, the effective-medium model 
is really accurate and reliable to use. In the present work, 
we again choose lithium tantalate (LiTaOH as the mate¬ 
rial of the wires like the previous study [23|. 

The dielectric function of LiTaOa is described by the 
Drude-Lorentz model H3- 


Sr — Sri 


jujj 


( 12 ) 


in which ujt/2'k = 26.7 THz, uji^/2tt = 46.9 THz, and 
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£oo = 13.4 |2^. This frequency dependence is in a 

very good agreement with the experiment. Thirty nine 
(39) THz is the frequency where the real part of the rel¬ 
ative dielectric constant is equal to 5ft(er) = —11.4, that 
is the target at the present stage. 

Next, Eq. [10] determines the condition for the second 
type of topological transition. Since fv<^l the pole and 
zero of the perpendicular permittivity are located closely 
on the frequency axis. Therefore P fa Z fa —1 and, 
consequently, Ej. fa —1. This value corresponds to the 
frequency range 45.5-46.0 THz. Thirty nine (39) THz 
and 46 THz are rather distant frequencies, and the cor¬ 
responding resonances of the Purcell factor do not over¬ 
lap. However, on the present stage it is worth to show 
that the radiation of the parallel dipole really experiences 
both these resonances. 

In the following, we report the numerical results for 
the Purcell factor of a finite sample of optimized WM for 
subwavelength electric dipoles either parallel or perpen¬ 
dicular. We performed full-wave simulations using the 
CST Microwave Studio simulator. For reliability we ap¬ 
plied two solvers (FDTD and frequency domain ones) and 
attained the coincidence. In our simulations, the sample 
of the WM is always of cubic shape with same internal 
geometry though its dimensions were varied. A subwave¬ 
length electric dipole (the same as in our work [2^) is po¬ 
sitioned symmetrically - centered within the gap between 
four central nanowires. This location is nearly equivalent 
to the averaging over different positions of the dipole . 
The simulated structure is shown in Fig. |4l 

We calculated two types of the Purcell factor. The first 
one is the usual (or full) Purcell factor i.e. an overall en¬ 
hancement of the radiated power of the given dipole mo¬ 
ment in the environment different from free space. The 
second factor is the so-called radiative Purcell factor that 
is, indeed, the enhancement of the power radiated by 
the given dipole into free space. Definitely, for an opti¬ 
cally large sample of lossy WM the second factor is much 
smaller because the most part of radiation is absorbed in 
the medium. We prove by extended simulations that the 
obtained values for the full Purcell factor can be identi¬ 
fied with those of the Purcell factor of the infinite WM. 
As to the radiative Purcell factor, it is not so, and the 
results depend on the sample size. 


A. Parallel dipole emitter 

Geometry and orientation of the dipole is shown in 
Fig. 01 The full Purcell factor versus frequency is pre¬ 
sented in Fig. [5] The maxima of the Purcell factor have 
nothing to do with the finite size of the sample and can 
be identified with the Purcell factor of an inifnite WM. 
We checked the absence of dimensional resonances simu¬ 
lating the sample with three different sizes corresponding 
to 14 X 14, 16 X 16 and 18 x 18 nanowires. Since we keep 
the cubic shape of the sample, the length of the wires 
was in these simulations equal to 2.664 ^m, to 3.064 ^m 



FIG. 4. Sub-wavelength parallel electric dipole in the cubic- 
shape sample of wire medium. The color spot illustrates the 
radiation which attenuates mainly due to optical losses. 



Frequency [THz] 


FIG. 5. Full Purcell factor versus frequency for different sizes 
of wire-medium sample. Here, the radius of each wire and 
the lattice constant of the medium are 32 nm and 200 nm, 
respectively. The emitter is parallel to the optical axis. 

and to 3.464 /rm, respectively. Two peaks of the full Pur¬ 
cell factor located at 38 THz and 46 THz (i.e. very close 
to values predicted by our theory) are kept for all three 
sizes of the sample. Also, the values of the Purcell fac¬ 
tor at these maxima was kept. At the first topological 
resonance, the full Purcell factor is nearly equal to 500, 
and at the second one it is equal to 180. This stabil¬ 
ity means that dimensional (Fabry-Perot) resonances at 
the frequencies of our interest are absent. Really, calcu¬ 
lations of the optical decay using the effective-medium 
model shows that these resonances are damped by losses 
even for a structure with 12 x 12 nanowires. 

As to the second topological transition, it is important 
to analyze the resonance of the perpendicular permittiv¬ 
ity, which is different from that analyzed in the previous 
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(b) 


FIG. 6. a)~Perpendicular component of permittivity with re¬ 
spect to the frequency for lithium tantalate wire medium, b)- 
Parallel component of permittivity with respect to the nor¬ 
malized wave vector at / = 45.9 THz. The solid lines show 
the real part and the dashed lines determine the imaginary 
part. The blue, magenta, black and red colors correspond to 
the radius = a = 32 nm, radius = 1.25a, radius = 1.5a and 
radius = 1.6a, respectively. 



(a) 



FIG. 7. Full Purcell factor (FPf) with respect to the fre¬ 
quency for several values of the wires’ radii. The lattice con¬ 
stant is 200 nm, and the initial radius is a = 32 nm. The 
emitter is parallel to the optical axis and the sample includes 
18 X 18 nanowires. 


subsection due to much higher optical losses in lithium 
tantalate. The frequency dispersion of e± is shown in 
Fig. 6(a) The blue curve corresponds to the initial value 
a =32 nm and to Fig. [S] Due to high losses, the real part 
of ex does not change the sign at the resonance. At 46 
THz 5ft(ex) is minimal and smaller than unity. This is 
sufhcient to mimic the regime ex = 0 which offers the 
second topological transition. As shown in Fig. |6(a)[ in¬ 
creasing the radius of the wires (the lattice constant is 
fixed) causes a stronger resonance which makes 5R(ex) 
become closer to zero. If 5i(ex) crosses the zero axis, 
then the zeros of 5ft(ex) occur at frequencies which are 
close to one another and in this range the WM is not very 
lossy. Really, in Fig. 6(a) one may see that the maxi¬ 
mum of the lossy factor |55(ex)| occurs at a frequency 


slightly lower than the band in which 5i(ex) is minimal. 
Therefore the increase of a is favorable for the second 
resonance of the Purcell factor. Indeed, increasing the 
fraction of nanowires we risk to lose the first topological 
transition. Therefore, for given b there must be optimal 
values for the wire radius. This optimum is to be found. 

In accordance to the insight that the Purcell factor 
of a very anisotropic medium is determined by the ratio 
e||/£-Li we inspected also the parallel component of the 
effective permittivity. Figure. [6(b)| shows e|| versus nor¬ 
malized parallel wave number at 45.9 THz. We see that 
at this frequency, 5i(£||) is approximately uniform versus 
/3, and it is a positive value in between 0.6 and 0.85 for 
wire radius varying within 60% range. The increase of 
the wire radius is important for the perpendicular per- 
































mittivity where it allows the 2nd topological transition 
when the wire radius becomes larger than 1.6a. For the 
parallel permittivity in this frequency range the wire ra¬ 
dius is not very important. 

The optimization of the WM implies that we need to 
approach two resonances of the Purcell factor and to 
enhance the second resonance compared to the values 
achieved in Fig. [S] Increasing the radius of nanowires 
we pass through its optimal value. In Fig. 7(a) we show 
the simulated Purcell factor for wire radii increasing from 
a = 32 nm to 1.6a = 51.2 nm. With these values of the 
wire radius we only mimic the second topological transi¬ 
tion. However, we keep the 1st one. The Purcell factor 
corresponding to this first resonance in Fig. |7(a)| keeps 
nearly the same for different wire radii. However, the 
resonance frequency moves versus the wire radius and 
approaches to 46 THz. Two resonances start to overlap 
when the radius exceeds a. The best case corresponds to 
the wire radius 1.6a = 51.2 nm when FPp > 300 at 40-46 
THz. The two resonances can be still distinguished but 
are almost identical and the dip between them is small. If 
the target is the broadest frequency band in which Pur¬ 
cell factor exceeds 300, the best result corresponds in this 
plot to 1.6a. 

However, the wire radius 51.2 nm is also the lower 
bound of radii for which 5i(e_L) changes the sign, and 
allows the second-type topological transition to really 
occur. Therefore, the Purcell factor at this reson ance 
grows versus wire radii as it is shown in Fig. 7(b) As 
seen in this figure, the first resonance becomes sacrificed 
if the wire radius exceeds 1.75a. On the other hand, Pur¬ 
cell factor at the second resonance still increases versus 
the radius and the best result corresponds to the radius 
2.5a =80 nm. In this case the maximal Purcell factor is 
equal 2 • 10^ and this value is by an order of magnitude 
higher than the values obtained for WM in (2^ and for 
stacked media in [^ . This amazing radiation enhance¬ 
ment is revealed for a parallel dipole which was previously 
considered as a very inefficient emitter for WM. 


It is worth to notice that at the second resonance 
the ordinary (TE) mode may contribute into the Pur¬ 
cell factor. For small fraction volumes and at lower fre¬ 
quencies the perpendicular component of the effective 
permittivity is approximately equal to unity implying 
Pp (TE) Ri 1. However, at the frequency of the second 
resonance 5J(e_L) ~ 0 or even 5?(e_L) = 0. Therefore, at 
this frequency the TE mode may also play a noticeable 
role in the Purcell effect. In this work, we, however, do 
not share out the contribution of the TE-modes into Pur¬ 
cell factor and leave this question for our next paper. 

Our next plot - Fig. |8(a)| - shows the radiative Pur¬ 
cell factor. It is still plotted versus frequency for different 
wire radii. Comparing Figs. |7(a)1 and |7(b)] to Fig. |8(a)[ we 
see that the radiative Purcell factor feels the first topo¬ 
logical transition (although keeps very small compared to 
the full one). However, the frequency where RPp is maxi¬ 
mal is determined in a contest of the resonance for the full 
Purcell factor (total radiated power) and the requirement 



32 34 36 38 40 42 44 46 48 

Frequency [THz] 

(a) 



(b) 


FIG. 8. Radiative Purcell factor {RPf) with respect to the 
frequency. (a)-i?PF for several values of the wires’ radii. The 
lattice constant is 200 nm, and the initial radius is a = 32 
nm. The emitter is parallel to the optical axis and the sample 
includes 18 x 18 nanowires. {h)-RPF for different sizes of the 
sample. Here, the lattice constant is as before and the wire 
radius is 32 nm. 


of the impedance matching at the WM interfaces. This is 
why the frequency of this maximum varies when the wire 
radius increases (and, respectively, the wave impedance 
of the WM changes). 

Notice that around the frequency of the second topo¬ 
logical transition the radiative Purcell factor is smaller 
than unity due to the strong impedance mismatch (total 
internal reflection). Though the total radiated power is 
maximal in this frequency band the radiation into free 
space is lower than that of a dipole in absence of the 
WM. This is the result of the almost total confinement 
of radiation in the sample whose effective permittivity is 
smaller than unity (see Figs. |6(a)] and |6(b)[ ). 

Figure. |8(b)| clearly shows the absence of dimensional 
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resonances around the frequency of the first transition 
representing radiative Purcell factors for different sizes 
of the sample. Approximately all the peaks happen at 
the same frequency and peak value decreases versus the 
sample size in accordance to growing attenuation across 
the sample. In general, all our plots show that the full 
Purcell factor simulated for a finite sample can be really 
identified with that of the infinite WM. On the contrary, 
the radiative Purcell factor is determined by the losses 
in the sample and by the mismatch at the interface. At 
the first resonance the matching is quite good and the 
size reduction increases the resonant value of the radia¬ 
tive Purcell factor. However, the size of the sample can 
not be very small in order to keep the properties of WM. 
Practically, it implies that the optimal sample size cor¬ 
responds to 12 X 12 nanowires. Still, in this best case 
99% of radiated power are dissipated in the WM, and 
the radiative Purcell factor attains only 8.3. 


of polaritonic WM. 



(a) 


B. Perpendicular dipole 


Here, we report the results obtained for a dipole per¬ 
pendicular to the optical axis and compare these results 


with those for a parallel dipole. Figure 9(a) shows the 


full Purcell factor as a function of frequency for differ¬ 
ent radii of the wires. The first topological transition 
does not have impact to the radiation of the perpendicu¬ 
lar dipole. This result is in line with our speculations in 
[^ . The Purcell factor remains approximately uniform 
around the frequency of the first transition and its value 
in the order of magnitude matches the predictions of Q . 
Between 45 THz and 46 THz the Purcell factor has a 
local maximum which corresponds to the second type of 
the transition. As seen, the Purcell factor at resonant 
frequency increases as the wire radius increases. 

To see clearer the difference between the parallel and 
perpendicular dipoles, we have shown the full Purcell fac¬ 
tor for these two dipole orientations in Fig. |9(b)] when the 
wire radius is fixed as 1.5a =48 nm. At the first tran¬ 
sition the Purcell factor for the parallel dipole is higher. 
At the second transition the perpendicular dipole radi¬ 
ates more. For values of the wire radius within a and 2a 
there is the band in which the parallel and perpendic¬ 
ular dipole radiate approximately the same power. For 
the radius 1.5a it is the band 40-44 THz. Emitters will 
all orientations of their polarization radiate in this band 
nearly the same power. This is the manifestation of the 
nearly isotropic Purcell factor which is a surprising effect 
for a very anisotropic medium. This isotropy and the 
wide band of the Purcell factor allows us to suggest an¬ 
other application of the Purcell effect, namely the radia¬ 
tive cooling of microscopic objects. Thermal radiation is 
usually isotropic and broadband. Very high values of the 
broadband and isotropic Purcell factor imply the very 
efficient radiation of the excessive heat into WM, that 
may be important for micro- and nanolasers and other 
active devices which can be incorporated into a sample 



Frequency [THz] 

(b) 

FIG. 9. (a)-Full Purcell factor for several values of the wires’ 
radii. The emitter is perpendicular to the optical axis. Here, 
a = 32 nm and b = 200 nm. The sample includes 18 x 18 
nanowires. (b)-Full Purcell factor for two different orienta¬ 
tions of the emitter in the case of radius = 1.5a. 


IV. CONCLUSIONS 

We have theoretically studied two different topologi¬ 
cal transitions of the dispersion that may happen in the 
same wire medium at different frequencies. Previously, 
we already revealed one of these transitions for a special 
case when the WM is performed of polaritonic nanowires. 
In this work we have found the second type of topologi¬ 
cal transition which appears at higher frequencies where 
the perpendicular to the optical axis component of the 
effective permittivity tensor crosses zero. This regime 
can be also mimicked by the resonance of the perpen¬ 
dicular permittivity when its real part attains a positive 
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value much smaller than unity. Both these topological 
transitions correspond to resonances of the Purcell fac¬ 
tor for a parallel dipole. We have shown that optimizing 
the WM, one may bring the two resonances to one an¬ 
other. Overlapping the corresponding resonance bands 
we achieve the huge and broadband enhancement of the 
power radiated by a parallel dipole. We confirmed our 
theoretical expectations by extensive simulations. The 
two types of resonances for the Purcell factor correspond¬ 
ing to two topological transitions were specified for a sub¬ 
wavelength dipole located at the center of a finite-size 
sample of a polaritonic WM. We have studied the radia¬ 
tive Purcell factor and explained its significant difference 
from the full one. Finally, we compared and discussed 
the frequency behavior of the Purcell factor for dipoles 
oriented in parallel to the optical axis and perpendicu¬ 
larly. We have found the regime of the resonant isotropic 
Purcell factor. The physical effects theoretically revealed 
in this work allows the broadband and strong radiation 
enhancement which can be applied in optical nanosensing 
and in radiative cooling of small objects. 


We claimed all reported effects for WM performed of 
lithium tantalate. However, our effects can be also engi¬ 
neered with other polaritonic materials, e.g. crystalline 
or polycrystal silicon carbide (SiC). The work [s^ re¬ 
ported the fabrication of some polaritonic WM for the 
infrared operation grown in a dielectric matrix. Perhaps, 
there are no limitations for fabrication of polaritonic WM 
with free-standing parts (as in our simulations). How¬ 
ever, we believe that main results of the present work 
can be reproduced for nanowires in a dielectric matrix 
with host permittivity larger than unity. Free space as a 
host medium is chosen in this work only for simplicity of 
calculations. 
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